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We study in details the dynamics of the one dimensional symmetric trap model, via a real-space 
renormalization procedure which becomes exact in the limit of zero temperature. In this limit, the 
diffusion front in each sample consists in two delta peaks, which are completely out of equilibrium 
with each other. The statistics of the positions and weights of these delta peaks over the samples 
allows to obtain explicit results for all observables in the limit T — > 0. We first compute disorder 
averages of one-time observables, such as the diffusion front, the thermal width, the localization 
parameters, the two-particle correlation function, and the generating function of thermal cumulants 
of the position. We then study aging and sub-aging effects : our approach reproduces very simply the 
two different aging exponents and yields explicit forms for scaling functions of the various two-time 
correlations. We also extend the RSRG method to include systematic corrections to the previous 
zero temperature procedure via a series expansion in T. We then consider the generalized trap 
model with parameter a £ [0, 1] and obtain that the large scale effective model at low temperature 
does not depend on a in any dimension, so that the only observables sensitive to a are those that 
measure the 'local persistence', such as the probability to remain exactly in the same trap during 
a time interval. Finally, we extend our approach at a scaling level for the trap model in d = 2 and 
obtain the two relevant time scales for aging properties. 



I. INTRODUCTION 



A. Trap models as toy models to study aging phenomena 

Trap models provide a simple mechanism for aging |]| |2( . The basic phenomelogical idea is that the slow dynamics 
of glassy systems is governed by metastable states denned as "traps" in the coarse-grained configurational complicated 
landscape. The distribution of the energy of the traps is usually taken to be exponential 

p{E) = e{E)^e-^ (1) 

On one hand, this exponential form describes the lowest energies in the Random Energy Model |3| and the distribution 
of free energy of states in the replica theory of spin-glasses On the other hand, it appears for the largest barriers 
in the biased one-dimensional Sinai diffusion hj, |5| as well as in more complex disordered systems such as fractals 
and percolation clusters |(j , elastic manifolds U bubble dynamics in DNA [g] and sequence alignment algorithms . 
The ubiquity of this exponential form actually comes ^(J from the exponential tail of the Gumbel distribution which 
represents one universality classes of extreme- value statistics [11). 

The exponential density of energy corresponds for the Arrhenius trapping time r = e° E to the algebraic law 



with the temperature-dependent exponent 



«M = *(r > l)-fe (2) 



T 

1 L= 7f (3) 



At low temperatures T < T g , the mean trapping time J drrqir) is infinite and this directly leads to aging effects. This 
mechanism shows that the presence of broad distribution of trapping times (J2J is rather generic at low temperatures, 
since it simply emerges from the exponential tail of extreme- value statistics for the energy barriers. 



B. Previous results on aging properties in trap models 

The dynamics of the trap model has been studied in details in its the mean field versio n [T^. IT^. lb] } [pi as well 
as in the one-dimensional directed version, in relation with the biased Sinai diffusion pi flfi Il7l Il8l Il9( fecj. In both 
cases, aging properties are characterized by scaling functions of the ratio (t/t w ) of the two times involved. 
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More recently, it has been proposed in |2]J to study trap models on a hypercubic lattice in arbitrary dimension d 
with the following generalized dynamics : the particle can jump from site i to any of the 2d nearest-neighbor sites j 
with a hopping rate per unit time given by 

w i ^{a) = Yd e+ " aEi ~ K1 ~ a)Et (4) 
in terms of the parameter a S [0, 1]. The case a — represents the usual trap model where the rate 

depends only of the initial site via the trapping time Tj = e@ Ei distributed with J3J) : the particle spends at site i a 
time ti distributed with the exponential distribution of mean 

_!± 

fn(U) = ~e T * (6) 

n 

and then jumps with equal probability l/(2d) to one of its (2d) nearest-neighbor sites. Another interesting case is 
a = i where the rate Wi^j depend on the energy difference (Ej — E{). 

Monte Carlo simulations and scaling arguments [2l| have shown the possibility of a so-called "sub-aging" behavior 
for the probability + t w , t w ) of no jump during the interval [t w , t w + t] 

u(t + t w ,t w )^n(^\ (7) 

with an exponent v < 1 given in one dimension by |2l| 

v{a) = —— 8 
1 + [i 

This exponent was proven to be exact by mathematicians, first for the usual trap model a = |22|, and then for 
arbitrary a |2^|. On the other hand, the correlation function C(t + t w , t w ), defined as the probability to be at (t + t w ) 
in the same trap it was as time t w , was shown to present a "full aging" behavior 

C(t + t w ,t w )~C^(j-j (9) 

for the usual trap model a = [HI], and then for arbitrary a So there are two different time scales t w and t w 
which play a role in the aging of this model. Asymptotic forms have also been heuristically proposed and numerically 
tested in for n(< + t w ,t w ) and C(t + t w ,t w ). Finally, let us mention a recent interesting application : these 
properties of aging and subaging for the trap model are relevant to explain the numerical simulations on the dynamics 
of denaturation bubbles in random DNA sequences . 



C. Previous results on anomalous diffusion and localization properties 

Apart from aging properties discussed above, trap models are interesting for their anomalous diffusion and local- 
ization properties. In particular, in dimension d = 1, the averaged diffusion front is expected to take the following 
scaling form at large times [23, [2JJ |23 

¥ ^^^W) 9 "{W)) (10) 

where the characteristic length scale follows the sub-diffusive behavior 

(11) 



This exponent can be found via a simple scaling argument on Levy sums 25, 26] or by a real space block-RG analysis 
[2?| . However, the scaling function itself g(X) is not known, but has recently been studied numerically in [23| together 
with asymptotic behaviors proposed in the limit /1 — > 1. 
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Another important issue concerns the localization properties. The localization parameters 

-f oo 

Y k (t)= J2 ^Ml (12) 

n— — oo 

represent the disorder averaged probabilities that k independent particles starting at site at time in the same 
random environment are at the same site at time t. It has been proven in |22j that the limit 12(00) is strictly positive 
in the full domain < /i < 1. The values of the limits Ifc(oo) have been numerically studied in |23j with various 
approximations. 



D. Summary of main results 



The aim of this article is to provide a probabilistic description, sample by sample, of the symmetric trap model in 
the limit of very low temperature \i — ► 0. We have previously developed a similar analysis for the directed version of 
the trap model |2£|. Here in the undirected version, each site may be visited many times and this leads to essential 
changes. In particular, at time t, the important traps are the traps having a trapping time r > R(t), where the scale 
R(t) is not linear in t as in the directed version |2(ij|. but is sub-linear in time 



where T (p) may be expanded in fi as 



"«"-'isb> 



T„(„) » l-'-"[l+OH (14) 



The corresponding mean distance between these important traps is then given by 

m = &(ao (15) 

where the exponent agrees with previous studies (|11J) described above, and where the prefactor reads 

£ (/i) = l + O(/i) (16) 
In terms of these scales, we obtain the following explicit results in the limit /i — > : 
• Scaling function of the disorder averaged diffusion front l|10|l 

g,(X) = e-W / due-"— — + 0( M ) (17) 
Jn \X\+u 



Localization parameters (|12|) 

V,,.(n\ = 

(k+1) 

• Generating function of thermal cumulants 



Y k (p) = jr^r + O(fi) (18) 



Z„(s) = In < e~ s «Fi > = J dXe- x X ( — coth — - lj + 0{p) (19) 

The series expansion in s yields the disorder averages of rescaled thermal cumulants : the first ones are the 
thermal width 



c 2 ( M ) = t Hm ^ ^ = 1 + 0(/i) (20) 



and the fourth cumulant 



< n 4 > -4 < n 3 >< n > -3 < n 2 > 2 +12 < n 2 >< n > 2 -6 < n > 4 ^ . , . . 

c 4 (/i) = hm _ = -4 + 0(/i) (21) 
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and more generally 

c fc M = (2fc + l)!B 2fe + 0(/x) (22) 

in terms of the Bernoulli numbers B n . 
• Two-particle correlation function 



+00 +00 / I \ 

C{l,t) = V Vp(n,t|0,0)P(m J i|0,0)« Ii |„_ m | s r a (p)*,, + I ) (23) 

n=0m=0 ^ 

where the weight of the delta peak at the origin corresponds as it should to the localization parameter Y2 = 
2/3 + 0(/x) JTSJ, whereas the second part presents a scaling form of the variable A = |4r. The scaling function 
C M reads 

C M (A)= e - A ^ + 0( M ) (24) 



The probability Tl(t + f„,) of no jump during the interval [t w , t w + t] takes the scaling form Q 

R(t w ) 



ni/ + / »- ) = n,, [ // = =u„^g= [f (n)]^-^zj (25) 



with the scaling function 

= f 1 dz^e-'" (26) 
Jo 

In particular, we obtain the asymptotic behavior 

U(t + t w ,t w ) ~ f^srl [M + 0(Ai 2 )] (27) 

— I f+OO \ .1+^1 / 

• The probability C(t + t w ,t w ) to be at time (t + t w ) in the same trap at it was at time t w takes the scaling form 



O / + /,, 1 = C„ ( /, = Rl+ t {tw) ) =C»(h = foMj- ) I2S) 



with the scaling function which reads at lowest order in /i 



= <5„(fc) = / ^dzz^KKz) (29) 



In particular, we obtain the asymptotic behavior 

C(t + t w ,t w ) ~ ' [m + 0( M 2 )] (30) 

■j^-^OO \t w J 

We also extend the RSRG method to include systematic corrections to the zero temperature procedure via a series 
expansion in [i : the corrections of order fi of the observables described above are given in Appendix El We also 
extend our analysis to the generalized trap models (gj and obtain that the only observable that depends on a is the 
two time correlation U(t, f'), which takes the scaling form 

ni ^ t+t ^= n %( v =w^) (31) 

in terms of the result l|26|l for the a = case, but that otherwise all other observables described above are exactly 
the same as in the case a = 0. The reason is that the influence of a is purely local around a renormalized trap and 
does not change the renormalized effective model at large scales. Finally, we also extend our RSRG approach to the 
trap model in dimension d — 2 at a scaling level. 
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E. Organization of the paper 



The paper is organized as follows : 

In Section [HI we denned the renormalized landscape for the usual trap model a = and study its properties : 
in particular, we obtain the relevant length scale and the two relevant time scales. In Section ITTT1 we describe the 
effective dynamics in the limit [i — ► and compute one-time and two-times observables in this limit. In Section llVl 
we study the corrections to the effective dynamics at first order in /i and we describe the hierarchical structure of 
the important traps that play a role at order fi n . In Section we extend our approach to the generalized trap with 
parameter a E [0, 1]. Section IVll we extend our RSRG approach to the trap model in d = 2 at a scaling level. The 
conclusions are given in Section IVIII and the Appendices contain more technical details. 



II. DEFINITION AND PROPERTIES OF THE RENORMALIZED LANDSCAPE 



A. Notion of renormalized landscape at a scale R 



We wish to adapt the Real Space Rcnormalization procedure already defined for the Sinai model |5j] and for the 
directed trap model [2(j to the undirected trap model. The basic idea is that the dynamics at large time is dominated 
by the statistical properties of the large trapping times. The renormalized landscape at scale R is defined as follows : 
all traps with trapping time t% < R are decimated and replaced by a "flat landscape" , whereas all traps with waiting 
time Tj > R remain unchanged. The distribution of the distance I between two traps of the renormalized landscape 
at scale R reads 



Pr(1) = 



I - / drq{T) 
R 



l-l 



drqir) 



(32) 



where the first part [..] represents the probability that (I — 1) traps have a trapping time r, < R, and where the 
last part represents the probability that the I trap has a trapping time Ti > R. So the appropriate rescaled length 
variable at large scale R is 



x -rT< 

and the scaling distribution is simply exponential 

V(X) =e- A 

The distribution of the trapping times of the traps in the renormalized landscape at scale R is simply 

9(r) 



Qr(t) = 0(t > R) x 



T \ T 



(33) 



(34) 



(35) 



In the directed version of the model, the particle visits each site only once, and the RSRG analysis directly deals 
with the trapping times Tj. However here in the undirected version of the model, each site may be visited many times, 
and thus it is necessary to introduce the notion of 'escape time' as we now explain. 



B. Notion of 'escape time' from a renormalized trap to another renormalized trap 

We now introduce the notion of the 'escape time' T from a trap To existing in the renormalized landscape at scale 
R. This trap is surrounded by two renormalized traps that are at distances 1+ and /_ on each side (see Figure ^1 . 
Whenever the particle escapes from the trap r , it can escape on either side with probability (1/2). If it escapes on 
the left, it will succeed to reach the trap r_ with probability and if it escapes on the right, it will succeed to 

reach the trap t + with probability 1/Z+. Otherwise, it will be re- absorbed again by the trap tq. 
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FIG. 1: Definition of the escape time from a trap in the renormalized landscape : the trap of escape time to existing in the 
renormalized landscape at scale R is surrounded by two renormalized traps that are at distances 1+ and £_ on each side. The 
escape time To is the mean time needed to reach either r+ or r_ when starting at tq. 



1. Number of sojourns in a renormalized trap before escape to a neighbor renormalized 



As a consequence, the probability En(n) to need (1 + n) successive sojourns in the trap tq before the particle 
succeeds to escape either to the trap t_ or to the trap r+ reads 



E R (n) 



1 



1/1 1 



-I— - 

2 I + + L 



2 \l + l- 

For large R, since we have l± = R^X± l|33[) . the number n follows the same scaling : the rescaled variable 



w 



R^ 



is distributed exponentially 



where the coefficient 



£(w) =a(A + ,A_)e- toa ( A +< A -) 



a(A+, A_) = 



1 / 1 



2 V \, + A 



(36) 

(37) 
(38) 
(39) 



depends on the two rescaled distances to the next traps. Its distribution k(a) over the samples has for Laplace 
transform in terms of the Bessel function K\ (IE 111 



k(s) 



+ OC 



dae- sa k(a) 



+ 00 



dXP{\)e~^ 



2sK%(V2s) 



In particular, the non-analytic behavior at small s (|E. 



k(s) ~ l-sln-+0(s) 



corresponds to the following algebraic decay at large a 



k(a) ~ - 



(40) 



(41) 



(42) 
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The decay at large s (IE2 



k(s) ~ .^/-e-^' (43) 



corresponds to the following essential singularity at small a 



k( a ) ~ (44) 



2. Total time spent inside a renormalized trap before escape to a neighbor renormalized trap 

Let us now consider the probability distribution Pi n (Un) of the total time U n spent inside the trap tq before its 
escape. It can be decomposed into the number n of sojourns, where n is distributed with i|36[l 

1+71 



tin =J2 li ( 45 ) 



i=l 



where t t is the time spent during the sojourn i in the trap To, so it is distributed with the exponential distribution 
© with Tj = To. Actually, since n is large in the large R limit, we have the central-limit theorem 

t in — n <U>= nr (46) 

n — >oo 

This explains the numerical observation |23| that the results are unchanged if the particle spends a time exactly equal 
to Ti at each visit to site i, instead of a random time ti distributed with J^J. 

Since the number n is distributed with (|37l I38[) . we finally obtain that U n is also exponentially distributed 

tin 

PUUn)^ To (47) 



with the characteristic time 



To = -R^T (48) 



Since the smallest trapping times existing in the renormalized landscape at scale R is tq = R, the time spent inside 
the trap To before it succeeds to escape scales as 

Un - R 1+f * (49) 

R — >oo 

3. Total time spent during the unsuccessful excursions before the escape 

Among the n unsuccessful excursions, there are m excursions on the left and (n — m) excursions on the left, where m 
is distributed with the binomial distribution 2 _ ™C™. Since n and m are large, we again have a central-limit theorem 

m n—rn 

tout = XX 7 +^2tt- m<t - > ~K n -m) <t+> (50) 
i=i j=i 

where < t± > represents the mean time needed to return to when starting at 1 without touching the point l± in a 
flat landscape. The asymptotic behavior (see Appendix lAl for more details) 

< i± > ~ l f (51) 

gives that the scale of t out at large R reads 

tout - nl - R 2 » (52) 
which is negligible with respect to ti n (|49|l for /i < 1. 
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4- Time spent during the successful excursion to escape 



We finally consider the diffusion time tdiff of the successful escape to the neighbor renormalized landscape. The 
free diffusion over a length I ~ i? M takes a time of order l 2 (see Appendix 1X1 for more details) and thus the scale of 
t di ff at large R 

tdiff ~ R 2 » (53) 
is the same as t out 152fl but is negligible with respect to ti n l|49l) for fi < 1. 



5. Conclusion 



So we obtain that the total time 

tesc = tin + tout + tdiff (54) 

needed to escape is actually simply given by the time ti n spent inside the trap tq . So the distribution of t esc is given 
by the exponential l|47|l with the escape time Jo- 
in conclusion, a trap of the renormalized landscape at scale R has a trapping time r distributed with (|35|l . but has 
an 'escape time' proportional to r 

T=-R»t (55) 
a 

with a factor R^ that explains the occurrence of two different time scales in this model, and with a prefactor a 
distributed with gUJ) 



C. Distribution of 'escape times' in the renormalized landscape 



The distribution of the escape time T in the renormalized landscape at scale R reads 

r+OO r+oc i 

Q R (T)= dTq R { T ) dak(a)S(T - -ff'r) (56) 

Jr Jo a 

It thus presents the scaling form 



where the scaling function reads 



In particular, for large T, there is the same algebraic decay with index (1 + /i) as for the distribution of trapping 
times 

M f ) T fLffc c ^ (59) 

where the constant c(/Lt) may be computed from the Laplace transform (|40|l using i|E4fl 



c(m) ee / + °° da^ = JL dss ^~k( s) = 2 ^^&±^ = 1 - M (l + lE - In 2) + 0( M 2 ) (60) 



For small T, we use the asymptotic behavior of k(a) at large a (|42ll 



1 lna+ (ln2- 1 -7 B ) „Tna 

fe ( a ) ~n~ + 3 +°( — ) ( 61 ) 
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to obtain 

q (f) ~ -J!— + -V—f(lnl+ln2- lE -l + -}—) +o(f 2 lnl) (62) 

Using l|40[) . the most convenient way to characterize the scaling function in closed form is by the following transform 
in terms of a variable w 

/ dTe-^Q ll [T)= fidvv^ 1 dak(a)e-^ va = dzz 1+2ll K 2 (z) (63) 

Jo ^ ' Jo Jo w M Jo 

D. Choice of the renormalization scale R as a function of time 

For small /i, the probability distribution (t^J is dominated by its long tail (|59p. and we may approximate it by 

^ [ Toip) 
f \ f 



Q»(f)^6(f>W)t(?m) (64) 



where the cut-off Tq chosen to preserve the normalization is determined by the coefficient for the long tail part (|59I6U|) 



T (m) = (cfc))* ^ 2e- x -™ 
(J,— >o 



2 



(65) 



For the unrescaled probability distribution Ij57(l . this corresponds to the cut-off 

T (m) = fl 1+ "T (/i) (66) 

It is thus convenient to associate at time t the renormalization scale R{t) such that 

To (^) = t (67) 

meaning that at time t, only traps with escape times T > t have been kept, whereas all traps with escape times T < t 
have been removed and replaced by a flat landscape. This leads to the explicit choice 

i 

= ^7T (68) 



.To(p) 

The corresponding mean distance between traps reads at this renormalization scale (|33|l 

m = [R(t)r = (69) 

with the prcfactor 

M = (c(n)r rh =1 + M(1 + 1e- In 2) + 0( M 2 ) (70) 

III. EFFECTIVE MODEL AT LARGE TIME IN THE LIMIT /x -> 
A. Effective rules for the dynamics 

The prescription for the dynamics is as follows : 

At time i, the particle starting at the origin O will be at time t either in the first trap Af+ of the renormalized 
landscape at scale R(t) on its right or in the first trap M_ of the renormalized landscape on its left. The weight Pm + 
of the trap M + is given by the probability to reach M + before A/_ for a particle performing a pure random walk, so 
it is simply given by the ratio of the distances from its starting point 

m^MM =^ = ^ PD 

P,M_„ +1 (M-|0) =^ = j J± I - ,72) 
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This rule for the effective dynamics is consistent upon iteration. Suppose there are three consecutive traps : the 
trap M_ is at a distance Z_ from the origin on the left, the trap M + is at a distance l + from the origin on the right, 
and the trap M ++ is at a distance I from the trap M + on the right. Suppose that the trap M + is decimated before 
the traps M_ and M ++ . The new weights for the traps M_ and M ++ become 

Pm_ = P[M_M + ](M-\0)+p [M _ M++] (M_\M + )p [M _ M+] (M + \0) = = P[m_m ++ ](M-\0) (73) 

Pm ++ = P[ M .M ++ ](M ++ \M + )p [M _M + ](M + \0) = - + *~ - - = P[M _ M++] (M ++ \0) (74) 

and thus the rules Q72JI for the occupancies of renormalized traps are consistent upon decimation of traps in the 
renormalized landscape. 



B. Diffusion front 



In this effective model, the diffusion front in a given sample thus reads 



1 



-V 



(0) 



where the scaling function reads 

7>(o) (x) 



X, 



X- 



8(X + X- 



X 



n \ 



X_ 



x, 



X- 



S(X - X+) 



(75) 



(76) 



in terms of the two rescaled distances X± between the origin and the nearest renormalized traps. Since the joint 
distribution of the two rescaled distances is completely factorized 

V{X + ,X_) = 9(X + )6{X-)e- x +- x - 



we obtain that the scaling function for the disorder averaged diffusion front p0(l reads at lowest order in /i 



9 { °HX) - 

In particular, its asymptotic behaviors read 

9 { "\X) 

9 {0) (X) 



dX 4 



+ oo 



dX-V{X + ,X-)V (0) (X) 



-\x\ 



due 



X 



-\x\ 



\x\- 



•oo \X\ 

1-1X1 



|X|-0 



In 



1 

W\ 



IE 



O |X| 2 ln- 



\X\ 



(77) 
(78) 

(79) 
(80) 



It is interesting to compare with the simple exponential front e x obtained in the directed version of the same trap 
model. At infinity, the front is reduced by the power p^y with respect to the exponential representing the distribution 
of the distance to a renormalized trap, because of the probability l/l to escape to this trap instead of being absorbed 
by a nearer renormalized trap on the other side. On the other hand, near the origin, the front is enhanced by the log 
term because it is more probable to be absorbed by a trap which happens to be very near instead of being absorbed 
by the renormalized trap on the other side. 



C. Localization parameters 

In a given sample, the localization parameters read at this order 

k - Pm+ 



yr=p k M + +p k M. (si) 



and thus averaging over the samples yields 

y (0) = / dX+ / dX_V(X+,X)— -t — r = r (82) 

k Jo Jo '{X + + X-) k (fc + 1) V ' 

The agreement with the numerical simulations of |2^| obtaining Y 2 — > 2/3 and Y 3 — > 1/2 is the clearest numerical 
evidence of the validity of the RSRG effective dynamics with only two relevant traps in each sample in the limit 
H -> 0. 
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D. Thermal width 

In a given sample, the thermal width reads at this order 

< An 2 (t) > ( °)= (n M+ + n„_) 2 pM + PM_ = n M+ n M - (83) 
and thus we obtain after averaging over the samples 



< An 2 (t) > r +oc f + °° 
c 2 {fi -» 0) = lim V = / dX+ / dX-V{X+, X_)X+X- = 1 

E. Disorder averages of thermal cumulants 

The generating function of disorder averages of rescaled thermal cumulants reads 



(84) 



Z{s) ee In < e~ s «w > (85) 

thus reads at lowest order 



f-oo r+oo / xr v 



/- +00 / 1 

Z {0 \s) = dX+ I dX-V(X+,X-)\n( 

Jo Jo V x + 







A_ X++X- J y ' 

dAe~ A A ( y coth y - 1 J (87) 

The series expansion in s gives the disorder averages of the rescaled thermal cumulants (|22|) . 

F. Two-particle correlation function 

The two-particle correlation function reads 

+oo +oo / , 

C ^t) = E E P(n)P{m)6 Qn _ m \ ^y^Sifi + -r-rC^ A = — 

n=0m=0 ' ^ ' ^ ^ ' 

where the weight of the delta peak corresponds as it should to the localization parameter Y^ ' = | discussed above, 
whereas the scaling function of the long-ranged part reads at lowest order 

C j? ) W=J dX+ ^ riA_P(A + ,A_) ry ±-^ Sx , x++x _ = e- x - (89) 



(A++A_) 2 — - 3 



G. Aging and Sub-aging properties 

As ex plai ned in the introduction, there are two different correlation functions which present different aging properties 
|22l l23l |24j|. We now very simply recover within our framework the expected sub-aging and aging exponents, and 
moreover compute the scaling functions in the limit T —* 0. 

1. Probability Il(t + t w ,t w ) of no jump during he interval [t w ,t w + t] 

The probability Tl(t + t w , t w ) of no jump during he interval [t Wl t w + t] is directly related to the probability ipt„ (j) 
to be at time t w in a trap of trapping time r via the integral 



n{t + t w ,t w ) = / 

Jo 



drip tw (T)e-i (90) 
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In our approach, we have assumed as a starting point that at lowest order in /it, the probability iptij) was given by 
the trapping time distribution in the renormalized landscape at scale R{t) I|tj8|) 

Wr)= te(f) (r),9(r>ii(i))^ (91) 



So we obtain the scaling form Q 



t \ ^ I r ~ , .. i t 

R(ty 



H(t + t w ,t w ) =IlJg = — — 4 j= [Tb(/x)]^— (92) 



and the scaling function reads at lowest order in 

flW(ff)= [ dznz»- 1 e-'° (93) 
In particular, its asymptotic behaviors reads 



n<%) ^ o i-. 9T ^ + (. 9 2 ) (94) 
nW(j) * *£±d (95) 

2. Probability C(t + t w ,t w ) to be at time (t + t w ) in the same trap at it was at time t w 

According to the analysis of the escape time from a renormalized trap to a next nearest one, the total escape time 
(|54|l is actually dominated by the time ti n spent inside the renormalized trap. As a consequence, we have 

C(t + t w ,t w ) = J dT<p tu (T)e-* (96) 

where <pt m (T) represents the probability to be in a trap of escape time T in the renormalized landscape at scale R(t w ). 
At lowest order in 4>t(T) is simply given by the probability Qrh)(T) given in l|57|l . So we obtain the scaling form 



(97) 



Cit + t w ,t w ) = <7„ (h = = C, (h = f (,)±) 

where the scaling function reads at lowest order in \x 

C^h) = / dTQ„(T)e-T (98) 
Jo 

in terms of the function Q M (T) introduced in (|58(l . Using 163(1 , we thus obtain the scaling function in terms of the 
Bessel function K\ as 

2u 

C ^ h) = j2^J dzz + flR i( z ) (") 
In particular, the asymptotic behaviors read 

^l-^Mni + OW (100) 

W - TJ ^A (WD 

IV. CORRECTIONS TO THE ZERO TEMPERATURE EFFECTIVE DYNAMICS 

In the previous Section, to compute all observables in the limit /i — > 0, we have considered that the distribution of 
the escape time was infinitely broad in the following sense : all traps with Tj < t were such that 5- ~ 0, whereas all 
traps with Tj > t were such that -j- ~ +oo. For finite /Lt, we have to take into account that these ratios are not really 
zero or infinite. This can be done in a systematic procedure order by order in \i as we now explain. 
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FIG. 2: Construction of the important traps in a given sample sample for a particle starting at the origin. The dashed line 
separates the "small" traps (that have a trapping time smaller than R(t)) from the "big" traps (that have a trapping time 
bigger than R(t)). The first big traps on each side called M+ and M- are occupied with a weight of order 0(/jP). The next 
big traps M++, M and the biggest small trap S in the interval ]M_, M+[ are occupied with weights of order 0(h). 

A. New effects at order y, 

At first order in /j,, we need to consider the following effects (see Figure [2} : 

• A trap M of the renormalized landscape has a escape time Tm which is not infinite. There is a small probability 

_ t 

(1 — e t m ) that the particle has already escaped from this trap at time t. 

If it has escaped, it has been absorbed by one of the two renormalized neighbors, with probabilities given by 
the ratios of the distances. We will say that the particle is "in advance" with respect to the effective dynamics 
of the limit fi — ► 0. 

• The biggest trap in an interval between two renormalized traps that we will call S , has an escape time T$ < t 

t_ 

which is not zero and thus there is a small probability of order e T s that the particle is still trapped in S at 
time t. We will say that the particle is "late" with respect to the effective dynamics of the limit fj, — > 0. 

The first corrections of order [i for various observables are given in Appendix El From a technical point of view, it 
turns out that the computation of averages over configurations becomes rapidly much more involved than the similar 
computations for the directed case [2(j ■ This is due to the fact that the escape time depends on both the trapping time 
and the distances to the neighbor traps. So here, contrary to the directed case where we have computed corrections 
up to order /i 2 pflj . we will not go further into the explicit corrections of higher orders, but simply describe the 
hierarchical structure of the important traps that appear at a given order in /i. 

B. Hierarchical structure of the important traps 

The procedure that we have described up to order \x can be generalized at an arbitrary order n as follows : all 
observables at order //" can be obtained by considering a dispersion of the thermal packet over at most (2 + n) traps, 
that have to be chosen among a certain number f2 n of possible configurations of the traps. Our aim here is simply to 
get some insight into the set of important traps that play a role at a given order n. 

At order /i™, the important traps are : 



• the main traps M_ and M + 
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• the next n large renormalized traps on each side : 
M 2+ ,...M( n+1)+ on the half line }M + , +oo[ 
M/ n +l)-,...M2- on the half line ] — oo, M_[ 

• the n biggest traps S^.-.Sn^ among the small traps in the interval ]Af_, M + [ 

• the (n — 1) biggest traps S^ .-Sn among the small traps in the interval ]M + , Af2+[, and the (n — 1) biggest 
traps S 1 ^ ...Sn among the small traps in the interval ]Af 2 _, A£_[ 

+ (2) +(2) 

• the (n — 2) biggest traps S 3 ...Sn among the small traps in the interval ]Af 2 +, M3+[, and the (n — 2) biggest 
traps S 3 ...S n among the small traps in the interval ]M 3 -, Af 2 _[ 

• ... 

• the biggest trap Sn^ n among the small traps in the interval ]M(„_x)+, Af n+ [, and the biggest trap S n ^ 
among the small traps in the interval ]Af„_, Mr n _n_[ 

The total number of traps is thus 

n 

T„ = 2 + n + 2^i = 2 + n(n + 2) (102) 

i=l 

which generalizes the number Tq = 2 (Af_,Af + ) at order n = and the number T\ = 5 (Af_,M + ,M jAf+^S 1 ) at 

order n = 1 

With these T ra traps, one has to construct the possible fl„ configurations of (2 + n) traps, that are ordered in 
positions, and that contribute up to order /x™, as explained in more details in Appendix iDl 

C. Discussion 

In conclusion, the effective dynamics for the trap model valid in the limit \x — > is also a good starting point to 
study the full aging phase < fi < 1, since one can build a systematic series expansion in /z for all observables. 
Moreover, the hierarchical structure of the traps that are important at order n give a clear insight into the structure 
of the dynamics for finite fi. 

V. PROPERTIES OF THE RENORMALIZED LANDSCAPE FOR THE GENERALIZED MODEL 

We now consider the generalized model defined by the hopping rates (@J with parameter a £ [0, 1] and study the 
changes with respect to the usual trap model corresponding to the special case a — that we have studied in the 
previous Sections. 

A. Definition of the renormalized landscape 

Here the renormalized landscape is defined for the energy variable : at scale T, all traps with energy E < T 
are decimated and replaced by a flat landscape, whereas all traps with energy E > V remain unchanged. So the 
distribution of the energy in the renormalized landscape simply reads 

Pr (E) = 9(E>r)-L e -^r (103) 

1 g 
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B-2 B-i A B +1 B +2 



FIG. 3: Definition of the two escape times : a deep trap A of the renormalized landscape has for nearest neighbors B±\ and 
for next nearest neighbors B±2. In the case a > 0, the trapping time te of A alone depends on a, whereas the trapping time 
t% of the interval made with the three consecutive traps B-iAB+i does not depend on a 

B. Trapping-time of a renormalized trap 

Let us consider a trap called A with energy E > T existing in the renormalized landscape : it is surrounded by traps 
(_Bi,£?2) on the right and (B_i,i3_2) on the left which have E ~ (see Figure|3J|. As a consequence, the hopping 
rates to escape from A are 



so we may associate to the trap A the trapping time 

T E = e /5(l-«) E 

The distribution of te in the renormalized landscape at scale T reads 

/j,(a) /r r \^( Q ) 



Tr ( T ) = °( T > T r)— — {—) 



with the cut-off 



and with the new exponent 



0(l-a)T 



l-a 



(104) 
(105) 
(106) 
(107) 
(108) 



instead of 



C. Trapping-time of the vicinity of a renormalized trap 



The hopping rates from the site B\ ( or i?_i) nearest neighbor of the renormalized trap A of energy E read 

1 

2 
1 

2 



(109) 
(110) 
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so there is a very high probability to return immediately to the deep trap A. 

We are interested into the probability distribution P* (t* ) of the total time t* spent at A before the first passage at 
£?2 or B_2 when starting at A. We may decompose t* into the number k of sojourns at A needed 

k 

** = 5> an) 

z=l 

where ii is the time spent in A in one sojourn, so it is distributed with f TE (ti) @- Denoting p2 the probability to 
hop to B 2 when starting at B\ 

_ W Bi ^b 2 _ 1 -f3aE n , 9) 

P2 w Bl ^ B2 + w Bl ^ A ~ 1 + e+^ E ~ 
we may express the probability distribution W(k) of k = 1, 2, ... as 

W(fc) = (l-p 2 ) fe -V 2 (113) 

The Laplace transform of t* (|lllfl thus reads 

f + OO +oo / 1 \ ^ 1 



/ dt*Pjt*)e- st ' =Y / W(k)( — 



st e J 1 + s% 



(114) 



so f* is exponentially distributed 



PJt*) = — e % (115) 

T E 

with the mean time 

r E = ^ = e fjE (l + e-^ E )^eP E (116) 

which is independent of the parameter a, and exactly coincides with the trapping time of the usual trap model (a = 0) 
distributed with @. In conclusion, as a grows from to 1, the mean sojourn-time t b <|1U5|) in the trap A decays, 
but this is exactly compensated by the large probability to return immediately back to A when on nearest neighbor 
sites B±\. 



D. Escape time from a renormalized trap 

Since we have found that the trapping time r* E (|116fl of the region B-iAB+i containing a renormalized trap A 
with its two neighbors was exactly the same as the trapping time t of the usual trap model (a = 0) studied in the 
preceding Sections, we obtain that the escape properties from a renormalized trap to another renormalized trap are 
exactly the same as for a — 0. So the distribution of the escape time T is given by JS7J, where the RG scale R(t) for 
the trapping times t* e is related to the energy RG scale T(t) via l|116|) 

R(t) = e mt) (117) 



E. Conclusions in d = 1 



All results based on the escape times T in the renormalized landscape do not depend on a : in particular, the 
diffusion front, the localization parameters, the disorder averages of thermal cumulants, and the correlation C(t' , t) do 
not depend on a. Among the observables we have considered, the only change will be for the probability TL(t w + 1, t w ) 
of no jump between t w and t w +t, since it directly involves the distribution of trapping time t b <j$0|l - At lowest order 
in /i we have that the probability distribution ip t (t) to be in a trap of trapping time r at time t is now given by l|106|) 

U r)= q £\r)=e [ r>r m )^(?f-y {a) (118) 
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where the cut-off is related to the RG scale R(t) JSBJl via p07H16fl 1 

r m = eK 1 -*™ = [Rit)] 1 -" (119) 

So we obtain the scaling form 

n( ^ t+t ^= n ^{ v =m^) (120) 

in terms of the scaling function i|93[l with the modified exponent I|1U8|) . 

F. Effect of a in arbitrary dimension d 

It is clear that the analysis presented above in d = 1 can easily be generalized for the trap model on a hypercubic 
lattice in arbitrary dimension d : a very deep trap A of energy E has a trapping time te (I105f) distributed with an 
algebraic law of exponent n{a) = /i/(l — a) <|1U8|) . whereas the cluster made out of A and its (2c?) nearest neighbor 
sites has for trapping time t e ~ e^ E <|116|) distributed with an algebraic law of exponent /i(0) = \x. This explains 
why the quenched model and the annealed model remain different even in arbitrary dimension d when a > 0, as was 
found in the numerical simulations |2lj | : in the quenched model, the exponent \x in the distribution of t* e comes from 
the very high probability to return immediately to a deep trap when leaving it, and thus in the annealed model where 
these returns are absent, the only relevant exponent is = ^/(l — a). As noted in |2l|, this effect is particularly 
clear in the region < /i < 1 < j^— : the annealed model is above its glass transition, whereas the quenched model 
is below its glass transition. 

In conclusion, in the glassy phase ji < 1 in arbitrary dimension d, the diffusion front, the localization parameters, 
the disorder averages of thermal cumulants, and the correlation C(t' , t) do not depend on a, and the only observables 
sensitive to a will be again those that measure the 'local persistence', such as the probability Yl(t w +t,t w ) to remain 
exactly on the same site between t w and t w + t. This conclusion is in agreement with the qualitative argument given 
in Appendix B 3 of pl| . 



G. Generalization to arbitrary jump rates satisfying detailed balance 

It has been recently argued in [2^| how the decay time for the function H(t w + t,t w ) should depend on the form of 
the jump rates satisfying detailed balance 

e pE >w(Ei to Ej) = e pE 'w{E l to Ej) (121) 

when they are not of the special form |0J. It is thus interesting to consider this question from the RSRG perspective. 

The arguments presented in details above for the special form Q of the rates may be straightforwardly generalized 
as follows. The trapping time associated to a trap A with energy E > T existing in the renormalized landscape (|1()5|) 
is more generally given by 

te = —r^- r (122) 

w{E to 0) V ' 

whereas the trapping time of the vicinity of this renormalized trap l|116|) reads more generally 

r * = te = w{0 to 0) + W (0 to E) 

Te p 2 w(E to 0)iu(0 to 0) { ' 

The detailed balance condition (|121|l for the special case Ei = E and Ej = thus gives 

t e (124) 



w(0 to 0) w(E to 0) w(0 to 0) 

So the condition to have exactly the same properties on large scales as in the usual trap model l@J with a = 0, is 
simply 

lim (e- 0E T E ) < +oo (125) 

E— »oo 
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since this condition is sufficient to get t* e ~ e@ E . In this case, the only observable sensitive to the precise form of the 
rates will be again the probability IL(t w +£, t w ) to remain exactly on the same site between t w and t w + t. Indeed, this 
probability will be a scaling function of the ratio {t / r w {t w )) , where T w (t w ) is the cut-off Tr(t m ) °f the renormalized 
landscape at time t w , i.e. at scale R(t w ) l|68|) . In terms of the energy cut-off T(t w ) 11171) 

\nR(t w ) \nt w 

T{tw) = —p~ ~ ^TT^) (126) 

we thus obtain 11221) the general prediction 

w(r(* w ) to o) w to 

Vice versa, if one is interested into constructing a model with a given decay time decay time [28| T w (t w ) for the 
function U(t w + t,t w ), it is sufficient to choose a form for the rates that satisfies the following asymptotic behavior 

w(E to 0) - — . (128) 

In particular, to obtain the ultra-slow logarithmic aging dynamics t m (t w ) = (hit w ) x , the jump rates have to behave 
as the power law w(E to 0) ~ E~ x for E — > oo, in agreement with |28j. 



VI. SCALING ANALYSIS FOR THE TRAP MODEL IN DIMENSION d = 2 



Since the random walk on a square lattice in dimension d — 2 is recurrent, we expect that here also, the scale of 
"escape times" in the renormalized landscape will be different from the scale of trapping times. On the contrary, for 
d > 2, the random walk is not recurrent, and there will be only one time scale given by the trappings times. 



A. Renormalized landscape 

As before, we define the renormalized landscape at scale R as follows : all traps with trapping time Tj < R are 
decimated and replaced by a "flat landscape" , whereas all traps with waiting time t, > R remain unchanged. 

We now consider the probability Pr(1) that the renormalized trap closest to the origin is at a distance / from the 
origin. At large scale, this means that the disc of surface S — irl 2 is empty of renormalized traps, whereas the disc of 
radius I + dl is not, yielding 

P R (l) ~ ^e~i?" (129) 

R^oo RV 

that generalizes I|33I34|) . Given a renormalized trap To, Pr(1) also describes the distribution of the distance I of the 
nearest neighbor renormalized trap. 



B. Escape time for a renormalized trap in d = 2 



We thus define the escape from a renormalized trap tq as the event where a particle starting at tq reaches for 
the first time the circle of distance I of the nearest renormalized trap : of course, contrary to the one-dimensional 
case, this does not mean that the particle has been really absorbed, but we expect that the particle has then a finite 
probability to be absorbed by a renormalized trap different from to. 

The probability to escape to the distance / without returning to To when starting nearby behaves as c/(lnZ), where 
c is a constant (see Appendix 151 for more details), and thus the number n of sojourns at To before the particle succeeds 
to escape scales as 

n~lnZ~lnJ?= (130) 
So the time spent inside the trap To before the escape scales as 

t m ~ nT - t In R% > R In R% (131) 
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The mean time < t >i of one unsuccessful excursion scales as (see Appendix [5] for more details) 

I 2 

<t>i ~ c 7 — 5- (132) 

so that the total time t out during the unsuccessful excursions scales before the escape scales as 

I 2 

tout ~ n<t >i~ — ~ - — —jr (133) 
In I In R.2 

It is thus negligible with respect to t in (|131fl at large scale. 

Finally, the time t^iff of the successful excursion over the length I scales as 

t diff ~l 2 ~R» (134) 

and is thus also negligible with respect to fj„ (|131|l at large scale. 

In conclusion, the total time needed to escape from a trap tq of the renormalized landscape at scale R 

tesc = tin + t ou t + trfiff {135) 

is dominated by ti n alone at large scale, and thus we may associate to the trap tq an escape time 

T = ar In R% (136) 

where a is a random variable of order 1 that characterizes the geometry of the neighboring renormalized traps around 
to- 

C. Choice of the renormalization scale 7? as a function of time 

The distribution of the escape time T 113(jfl in the renormalized landscape at scale R will thus present the scaling 
form 

%(r) = ^M f = 7^) <137) 

So the renormalization scale R has to be chosen as a function of time with the scaling 

t - R(t) In R% (t) (138) 

i.e. by inversion at leading order 

R(t) ~ -Lr (139) 

So the recurrence properties of the random walk in d = 2 induces again a different scaling for trapping times and 
escape times in the renormalized landscape. As a consequence, the correlation H(t + t w ,t w ) will be a scaling function 
of t/R(t w ) with (|139l) whereas the correlation C(t + t w ,t w ) will be a scaling function of t/t w , in agreement with the 
results obtained in |29ll30| . 

D. Number of important traps 

However, in contrast with the d = 1 case, there are two important length scales in d = 2 as we now explain. 
On one hand, given the choice H139|l for the renormalization scale R(t) as a function of time, the length scale for 
the distance between two renormalized traps is then given by i|129|) 



hit 2 



(140) 
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On the other hand, let us now recall the scaling analysis using Levy sums [2l|, |25|, |26| in d = 2 : after iV steps, 
the number of visited sites behaves as S = N/(lnN) and each site is typically visited In TV times, which yields the 
correspondence 



s 



=1 



N 



t~^nQnN)~QjiN)l — ) (141) 



Since the typical distance reached after N steps scales as r ~ v N, this yields in terms of the time t after the inversion 
of I141fl at leading order 



r(t) ~ n (142) 

So here, in contrast with the d = 1 case, the typical distance r(i) reached after time t does not coincide with the 
typical distance £(t) l|140|) between two renormalizcd traps. As a consequence, the number n(t) of important traps 
contained in the disk of radius r(t) will not remain finite but will grow logarithmically in time 

"(*) ~ jsS ~ h*** ( 143 ) 

This is in agreement with the absence of localization in d = 2 as t — > oo [3(j ■ 

However, in the regime where t — > oo and fi — > with i^ 1 fixed (|143|) . we expect that the diffusion front in a given 
sample may be described as a sum of delta peaks that are out of equilibrium with respect to each other, i.e. that their 
weights depend on the geometry of the renormalized traps in the region around the origin but not on their energies. 
It would be of course very interesting to build some RSRG effective dynamics for d = 2 case in this regime, but this 
goes beyond the scope of this paper. 



VII. CONCLUSION 



We have studied in details the properties of the one-dimensional trap model via a disorder-dependent real-space 
renormalization procedure. The RSRG approach provides the correct exponents for the relevant length scale and the 
two relevant time scales in the full domain < fi < 1, and then allows to compute scaling functions in a systematic 
series expansion in fi. Since we have already summarized our main results for disorder averages of observables in the 
Introduction II Dl we will not describe them again but instead discuss the physical meaning of our construction. 

We have seen that at lowest order in ^i, in a given sample, the particle can be only on two sites, which are the two 
nearest renormalized traps M + and M_ existing at RG scale R(t), and that the weights of these traps are simply 
given by the probabilities to reach one before the other one. This means that the dynamics remains out of equilibrium 
forever : the weights are not given by Boltzmann factors, they don't even depend on the energies of these two traps, 
but they only depend on the ratio of the distances to the origin! When including the first corrections in [i to this 
picture, we have taken into account the biggest trap S in the interval ]M-, M+[, and the next renormalized traps 

M and M ++ , but the dynamics still remains out-of equilibrium since the weights are still determined by the lengths 

between these traps and by escape rates of the form e~*/ T . This out-of-equilibrium character will thus persist in the 
whole localized phase < fi < 1 where these deep traps keep finite weights. This explains why the numerical studies 
pH |23| have found that the ideas of "partial equilibrium" and "effective temperature" were not able to capture the 
long-time properties of the trap model. This is similar to what happens in the directed trap model or in the related 
biased Sinai diffusion, where the thermal packet is also broken into sub-packets that remain out-of-equilibrium with 
each other even in the infinite time limit |20|: this is in contrast with the unbiased Sinai diffusion, where the thermal 
packet is asymptotically localized in one single infinite valley |3l| . in which particles are at equilibrium which each 
other [12. 

From the point of view of numerical simulations on random walks in random media, this shows that to study the 
localization and the convergence or not to a quasi-equilibrium regime at long times, it is interesting to study the 
dynamics in a single disordered sample, and not only disorder-averaged quantities : for instance, the pictures of a 
thermal packet in a given sample obtained in for the Sinai diffusion, and in [13 for the directed trap model, seems 
to give the clearer insight into the question of localization in one valley or in typically a few traps with finite weights. 
And whenever the structure of the thermal packet consists in a few sub-packets, whose positions and weights are 
sample-dependent, it is a very strong indication that some appropriate real-space RG description can be constructed 
to study the dynamics. 
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In a forthcoming paper |3J], we will show how the present RSRG approach for the unbiased trap model can be 
generalized to obtain explicit results in the limit fi — ► for the linear and non- linear response to an external bias, 
when it is applied from the very beginning at t = or after a waiting time t w . Recently, this problem has been studied 
via scaling arguments and numerical simulations in |35j , and was shown to satisfy a non-linear Fluctuation Theorem 
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APPENDIX A: STATISTICAL PROPERTIES OF EXCURSIONS IN d = 1 

As explained in the text, to study the excursions in the renormalized landscape (|5t)l53(l . we have to study the 
following standard problem : what is the probability distribution P x (t) of the time t of the first-passage at x = 
without having touched the other boundary x = I before, for a particle in a pure trap model? 

For x = 1, I — 1 this probability distribution satisfies the equation 

d t P x (t) = [P x +i(t) + P x -i(t) - 2P x (t)} (Al) 
with the boundary conditions Po(t) — 5(t) and Pi(t) = 0. So the Laplace transform with respect to t 

r+oa 

P x {s)= / dte- st P x (t) (A2) 
Jo 

satisfies 

P x +i(a) + Px-i(a) - (2 + s)P x {s) = (A3) 
for x — 1, I — 1 with the boundary conditions Pq{s) — 1 and P/(s) = The solution reads 



2 + s-Vs^+As 

PO) = g ^ ' 



In particular, the series expansion in s yields the first moments 



h(x)= f + dtt k P x {t) (A6) 
Jo 



For n — 0, the probability to reach before / when starting at x reads as expected 



For n= 1,2 



0o(x) = (A7) 



= x(l - x)(2l - x) 

h{x) _ x(l-x)(2l~x)(AP + 6xl-3x 2 + 5) (A9) 
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a. Unsuccessful excursions 

For the unsuccessful excursions present in t out (|50|l . we need to consider the special case x = 1, where the first 
moments read 

0o(l) = 1 - ~ (A10) 

W1) -tfflli! , ' (All) 

^(1) = C-l)(a-l)(4* 2 + 6Z-3 + 5) ^ il (A12) 
v ' 360Z z^oo45 v ; 

The leading terms can be understand as follows : the time t is of order I 2 with a probability l/l, which is the 
probability to arrive very near the forbidden extremity x = I. 

To get the leading behavior in I of all moments, we thus need to consider the rescaled Laplace transform 

1, r- , r- ~ ( 1 



Ag)^ = l- 7 (^coth^-l) + U) (A13) 



which corresponds to the following leading behavior at large I for the moments 

K (A14) 

6. Successful excursions 

For the successful excursion present in tdiff l|53(l . we need to consider the special case x = I — 1. Since in this case, 
the normalization is 

0o (J - 1) = 7 ( Al5 ) 



we consider the 'normalized' first moments 



~ ^ (A16) 



O (i - 1) i^oo 6 

(A17) 



9 2 (l - 1) . 7l 4 



9 (l - 1) z-oo 360 

To get the leading behavior in I of all moments, we thus need to consider the rescaled Laplace transform 

So the time tdiff for the successful diffusion scales as I 2 as expected. 

APPENDIX B: STATISTICAL PROPERTIES OF EXCURSIONS IN d = 2 

To study the scaling properties of excursions in the renormalized landscape in dimension d = 2, we have to study 
the following problem : what is the probability distribution P r (t) of the time t of the first-passage on the circle r = 1 
without having touched the other circle r = I before, for a particle starting at radius r and diffusing freely? 

For 1 < r < R, this probability distribution satisfies the diffusion equation in radial coordinates 

ftJV(t) = AP r (t) = + \ |) P r (t) (Bl) 



with the boundary conditions -Pi(t) = 5(£) and Pj(t) = 0. 
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So the Laplace transform with respect to t 

P r (s) = / dte- st P r (t) (B2) 
Jo 

satisfies 

^ + r*- S )^ ) = (B3) 

for 1 < r < R with the boundary conditions -Po(s) = 1 and P;(s) = 0. 

In particular, the moments appearing in the series expansion of P r (s) in s 

8 k (r) = / dtt k P r (t) (B4) 



may be computed by the recurrence 



\ dr 2 r dr 

For k = 0, the probability 6* (V) to reach first r = 1 before r = I when starting at r is the solution of 

d 2 1 d 



v dr 2 r dr 

with the boundary conditions 0q(1) = 1 and 8q(1) = 0, and thus reads 



> (r) = (B6) 



In t 

W = l-- (B7) 



For = 1, 0i (r) is the solution of 

d 2 1 d \ . , . / lnr 



rfr 2 ' r dr ) L ^ ^ \ hxl) ^ ^ 

with the boundary conditions 0i(l) = and 0i(£) = 0, and thus reads 

9l(r) = ^l,„ r+ ^-'){K- 1 ) (B9) 

4 In I 4 In / 

In particular, to study the unsuccessful excursions, we may consider r = 2 for starting point. The probability to 
escape is then 

*o(2) = 1 - ^ (BIO) 
and the mean time of an unsuccessful excursion has for leading behavior in I 

9l{2 \^ jf^l ln2 (BU) 
z^oo 4 in l 



APPENDIX C: EFFECTIVE MODEL AT FIRST ORDER IN fi 



In this Appendix, we compute the corrections at first order in /i by taking into account the two effects described in 
the text in Section ITVAI 
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1. Statistical properties for particles "in advance" at order 

There is a probability 

n a = (l-e-^) (CI) 

that a particle has already escape at time t from a renormalized trap of escape time Tm- 

We consider a particle starting at O that was in the trap M + in the effective dynamics of previous Section. With 

probability 7r a , the particle has already escaped from the trap M + at time t, and then it has been absorbed by M_ 

or by M ++ , the two nearest renormalized traps, with probabilities given by the ratios of the distances. 

We note as before X± the rescaled distances between the origin and M±, and we note A the rescaled distance 

between M + and M ++ . The joint probability distribution of X+, X,tm) is completely factorized 

Vm-M + m ++ {X-,X + ,\,t m ) = 0(X + )9(X_)0(\)9(t m > (^-Y e -x+-x—x (C2) 

Since the escape time Tm is proportional to the trapping time 

T M = —R^tm (C3) 

(IM 



with a prefactor om that depends on the positions via the prefactor 

1 / 1 

2 + A 
it is more convenient to use the variable 

A = X_ + X+ (C5) 

instead of X + . So the measure l|C2|l becomes 

dr M {R(t)Y*_ x 



V m _m + m ++ (X-, A , A, r M ) = 0(0 < X_ < X o )9(X)0(t m > R(t))»— e- A °- A (C6) 

To compute the average of an observable A(X-, Ao, A, tm) with respect to this measure, it is more convenient to 
first integrate over X^ € (0, Ao) and then to integrate over the remaining variables with the following notations 

A(X-,\o,\,t m ) =« [ °dX_A(X_, Ao, A, r M )» a (C7) 
Jo 

where the notation << .. >> a denotes 

r+oc j / r>n\ \ M ;-+oo r+oc 

«/(A ,A,r M ) »a= / dTM»— — / dA / G?Ae~ Ao ~ A /(Ao, A, tm) (C8) 

Jfl(t) r M V T M J Jo Jo 

For the simplest observables, we will need integrals of the following form, which can be computed in terms of Bcssel 
Functions via (IE1I) 



W M (/x)=«(l-e 3 ^U + ^i)AgA«» 



2^i 



r(l + p)r(l +q)- -^-K^^K^z) 



(C9) 



where the parameter zo depends only of /i via To introduced in the explicit choice of the renormalization scale R(t) 

(EH 



2(,(/i) = V T^W = ' (( lm 
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In particular, in the limit of vanishing //, Zo(ii) has for limit l|65|l 

z (0) = 2e~^ El = 0.90895 
so that the integrals read at lowest order in /i 



(Cll) 



W P M ~ 2/i 



20 (°) 



r(i + P )r(i + g ) - ^-^(z)^^) 



o( M 2 ) 



(C12) 



As first example, let us consider the probability to have already escaped from a renormalized trap at time t 



7r = 1 - e t a/ = 1 - e A o 
Its average over the configurations can be computed via 



TT a =<< A 



1 - e 



>>a 



(C13) 



(C14) 



which corresponds to the form (|C9|I for the special case (p = 1,5 = 0) and is thus of order fj, (|C12 



TTa = Wl,o(M) - 2 M 



(i+te) 

2e 3 dz 



0(/i 2 ) = /i 0.678238 + 0(/i 2 ) (C15) 
Let us now write the diffusion front in a given sample : it is a linear combination of the two possibilities I|72I74|) 

^M-m}m++ ( n ) = e"~ P[M-M+](M-\0)Sn,n M _ + P[M-M+] {M+ |0)£„,n M+ 



+(l-e t m) P[m_m ++ ](-W_|0)5„,„ m _ +P[m_m ++ ](M++|0)<5„,„ m++ 
and thus in rescaled distances, the correction with respect to the zero-th order (|76|l reads 



(C16) 



(0), 



(1 - e T « ) 



<5(X + XJ) + -^—5(X - (A + A - 



(C17) 



where the parameters are distributed with the measure (|C6|) . 

Similarly, if the particle has escaped from the M_ , it has been absorbed by M + or by M and the properties are 

the same as above by the symmetry X — > — X . 



2. Statistical properties of particles "in delay" 

We are now interested into the trap S with the biggest trapping time t$ in the interval ]M_, M+[. Its position is 
uniformly distributed in this interval, so the joint distribution of (X-, Xs, X + , ts) reads 



■D m _ sm+ (X-,X s ,X + ,ts) = 9(X + )9(X_)9(-X_ < X s < X+)9(t s < R(t))± fMV e -^+^-)(^) 



The measure for the positions alone reads 

J dr s V M _SM + (X-,X s ,X + ,rs) = 9(X + )9(X_)9(-X_ < X s < 
whereas the measure for the trapping time alone reads 

dX+ ( dX_ f dX s V M _sM + (X-,X s ,X + ,T S ) = 9{t s < R(t))2fx 



X- 



r 2/i-l 
J_s 

R 2 »(t) 



(C18) 



(C19) 



(C20) 
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Since the associated escape time 



T S = —R^ts 
as 



depends on the positions via the prefactor 



as 



2\(X + -X S ) (X_+X s ) 



(C21) 



(C22) 



it is more convenient to replace the rescaled distances (X_, X+) between M± and the origin by the rescaled distances 
between M± and the trap S 



x + = x + - x s 

X- = X- + X s 
With these new variables, the measure (|(J18|1 becomes 



V m _sm + ( x s,ts;\ + ,\-)=0(\+)6(\-)9(-\ + <X S < A_)0(r s < R(t)) — 



n (R{t)\»- {x++ x_)(m.y 



(C23) 
(C24) 



(C25) 



To compute the average of an observable A(Xs, A+, A_, ts) with respect to this measure, it is convenient to first 
integrate over Xs G [— A+, A_] and then average over the remaining variables 



A{X S ,X+, A-,r s ) =<< / dX s A(X s , A+,A_,r s ) » 

J — A i 



where the notation << .. >> d is defined by 



<< /(A+,A_,r s ) » d = 



R(t) 



j^(R(t)\ f „ f 



rs 



TS \ TS 



dA_ 



(C26) 



/(A+,A_, t s ) (C27) 



For the simplest observables, we will need integrals of the following form, which can be computed in terms of Bcsscl 
Functions via (|E1|) 

2, 



n Ptq (fi)=«e t s x p + x q _ » d = y 



where the parameter zq has been introduced in (|C10|I . At lowest order in \i we thus obtain 

2n '■+ 00 



/ dzz^+v+i) Kl+p {z)K 1+q (z) (C28) 

J z (/i) 



Let us first consider the probability 



MO) 



dzz^ + P+^K 1+p (z)K 1+q (z) + 0(ii 2 ) 



TT d = e T s = e 



' 2R( j -tq \ A 4_ A _ 



(C29) 



(C30) 



that a particle that has been trapped by S is still in the trap S at time t. Its average with respect to the measure 
(|C25ll can be computed as 

(C31) 



7^=«(A+ + A_)e »H%U+ + *y » d 
It is of the form (|C28|) and is thus of order [i (|C29|) 



n d = 2n x ,o(M) — 2/i 



H-oo 



2e 



dzz 2 K 2 {z)K 1 {z) + Q(p 2 ) ~ A1 1.20205... + O(m 2 ) 

2 ^ *^ 



(C32) 



Let us now write the diffusion front in a given sample at this order 

^m H sm + P0 = e"*0(X S > 0) [p[M_s](M_|0)<y„, nM _ +P[M_S](S|0)«„,„ fi 
+e _ *0(X s < 0) [p [SM+] (M + |0)5 n ,„ M+ +P[5M + ](<5|0)^„,„ S 

t r 

+(l-e T s) P[M_M + ](-^-|0)<5 n , nM _ +P[m_m+](-W+|0)5„,„ m 



(C33) 
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and thus in rescaled distances, the correction with respect to the zero-th order (|76|l reads 

= er^0(X s > 0) 



+e~^8(X s < 0) 



^S(X + (A_ - X S )) + {X - x Xs) S(X - X S ) 
(- Xs} 5(X (A + + X S )) + X ±±^6(X X S ) 



-e T s 



S(X + A_ - X S ) + X ~~ Xs 5(X - (A+ + X s j) 



A+ + A_ A+ + A 

where the parameters are distributed with the measure i|C25|) . 



(C34) 



3. Correction of order jj, to the averaged diffusion front 

a. Contribution of particles in advance 

To compute the contribution (|C17|1 to the diffusion front of the configurations M_M + M ++ , we first integrate over 
X_G[0,A o ] 



A 



dX-V 



(1) 

M_M + M + + 



(X) = (l-e ^7) 



Aq - \X\ 
A 



0{O<X < A ) 



MXl : 9(0 <-X<X )+ X + A ° 7 W 9(X < X < A + A ) 



A (A + A ) 



A + A 



(C35) 



So taking into account the configurations M Af_M + via the symmetry X — > —A", we finally get that averaging 

with the measure (|C6(1 yields the following contribution of particles in advance to the scaling function at order fi 



gW{X) =J dxj dX ^e-^) 
Jo Jo Jo w 



(l-e-^u+i)) 



_ 2 ^-^(0 < X < A ) + A + Ao ;l X l g(o < AT < A + Ao) 

Aq A + Aq 



(C36) 



b. Contribution of particles in delay 



To compute the average of the specific contribution l|C34|) of the configurations S-MS+ with respect to the measure 
2lj, we first integrate over X$ 



J dX s 8(-X + <X S < X-)V§_ SM+ (X) 



e J s 



e{\x\ < A_)^i + e { \x\ < x J x +7^ - .(ixi < a + + a ; A + A - A 



A_ M "X/ A 

so that the correction to the averaged diffusion front reads 

/>+oo f-\-oo /*+oo J 

9 < i ) {X) = J dX + dxj 
Jo Jo 

A- - |*| 



A+ + A_ 



(C37) 



e "^M^" 



*l < A_)^— ^ + 0(1*1 < A + ) (A+ : W) - < A + + A_) A - + A+ - A 'l 



A_ 



A^ 



A+ +A_ 



+ 0(^ 2 ) (C38) 
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c. Total correction at first order in [i 
The two contributions (|C36IC38|) yields the following total contribution at first order in fi 



/■+oo r+oo 

gW (X) ee 5 « (A) + gW (X) = fi d\ / dXe^^ 

Jo Jo 



IE + hi ■ 



Tb(0) ^ A + A 

An A 



-2- 



. An -|X| 



A 



0(0 < X < A 



A + Ap- \x\ 
A + A 



0(0 < X < A + A ) 



0( M 2 ) 



(C39) 



4. Correction of order /j, to the thermal width 

a. Contribution of particles in advance 
The specific contribution of the configurations M_M + M ++ to the thermal width reads with the measure l|C6|l 

(C40) 



After the integration over A_ l|C7|) . we thus obtain an integral of type l|C9|l for the values p = 2, q = 1 



(C41) 



By symmetry, the configurations M M—M+ give exactly the same contribution, and thus the total contribution 

of particles in advance reads at lowest order in \i (|C12|) 



[c 2 ]« = W 2 ,i(M) ^ W 



tl+7B) 

2e ~ 2 cfe 



4 -_A^)A 2 (z) 



0{n 2 ) = [i 0.5383 



(C42) 



b. Contribution of particles in delay 

The specific contribution of the configurations M_SM + to the thermal width reads with the measure (|C25|I 

^ [c 2 ]i?+ s ( ll + - [C2]^_ M+ = -e-*S [9{X S > 0)A+(A_ - X s ) + 6(X S < 0)A_(A + + X 5 )] 
After integration over As, we thus obtain the integral of type l|C28|) for the special value (p = 2, q = 1) 



rf 3 



;*~0 4 



<< e r s 

f + 00 



A 2 



+ A_- 



A 2 



>> d = -fi 2) i(/i) 



2 c 



(i±m 



2 2 

dzz 4 A 3 (z)A 2 (z) + 0(> 2 ) = -^2.2172... + 0(/i 2 ) 



(C43) 



(C44) 



c. Total correction of order for the thermal width 
Adding the contributions of particles in advance and in delay, we get the total correction at order 



6 



-/i 1.6789 



(C45) 
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5. Correction of order /i to the localization parameters 

a. Contribution of particles in advance 
For a given configuration M_M + M ++ , the localization parameters read (|C17p 



l(0) + (l) 



e j m 



X. 



i k 



(1-e t m ) 



i k 



A + A 



t_ X_ t_ X_ 

1-e t m- (l- e t m ) 

Ao Ao + A 



(C46) 



The average of the specific contribution may be computed with <|C7|) . so that taking into account the similar contri- 
bution of the configurations M M_M + , we finally obtain 



\ya (1) - 2lV i J (0)+(1) -fyj (0) "l 

[Ik\a — *\[ I k\M-M + M ++ V 1 k\M-M + ) 

Ao 



2 « 



-fe- 



1 + (1 - e t m ) 



k A o 



(A + A) fe 



Ed 



A + A c 



In particular, it reads for k = 2 



[r 2 ]« =« ^ 



2A 



(1 - e t m ) 



A 



(A + A ) 



-2 +(1 - e ^) 2 2 1 - 



AA 



(A + A ) 5 



(C47) 



(C48) 



b. Contribution of particles in delay 
For a given configuration M_SM + , the localization parameters read l|C34(l for the case X$ > 



l(0) + (l) 



k \M-SM + 



e T s 



A_ - X, 



A_ 



(1-e -s) 



A- - X s 
X+ + X- 



__^a_-a: s , _* a_-x s 

1 - e T s : - (1 - e " " 



X- 



\++X- 



(C49) 



and a similar expression for Xg < 0. 

The integration over X s <E [0, A_] l|C26j) yields 



A_ 



(fc + 1) 



e ^ + (1 - e"^) fc 



(A + +A_)^ 



m=0 



A+ +A_ 



(C50) 



The case Xg € [— A+, 0] leads to the symmetric contribution via the exchange between (A + , A_), and thus the average 
of the specific contribution reads 



[Yk]^ =« J_~ dX s ([Y k ] 



i(o)+(i) _ , v i(0) . 

\M_SM + l 1 k\M_M + I ^^d 



1 



= << 



(e" fe * - 1)(A_ + A + ) + (1 - e-^) k ' + ^ 



(fc + 1) 

In particular, it reads for k = 2 



(A++A_)* 



Ed 



(i) 



<< -e T s (1 - e r s ) 
3 v ; 



A_ +A, 



m— 1 



A_A_i 



A_ +A, 



e J s 



A_ A™ + A+AI 
(A+ + A_)'» 



(C51) 



(C52) 



6. Correction to Aging and Sub-aging properties 

a. Probability H(t w +t,t w ) of no jump during he interval [t w , t w + t] 

Since the probability Tl(t w + t, t w ) of no jump during he interval [t w , t w + t] is directly related to the probability 
if>t w (t) to be at time t w in a trap of trapping time r via (|90|l . the configurations M_M + M ++ do not give any correction, 
since the trapping time of M ++ has the same statistical properties as the trapping times of M_ and M + . 
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On the contrary, the configurations Af_5M + will give a correction to V't( r )- F° r the case X$ > 0, we have 
[Mr)mil + = e-^±^6(r - r s ) + (1 - ^^^(r - r M+ ) 



_^ A_ - X s A-Xs 

x ~ e s ^—- (1 ~ e s) a7Ta7 

so that the specific contribution yields after integration over Xs S [0, A_] 

A„ /-A 



5(r-r M _) (C53) 



dX 3 [Mr)]^_ SM+ = / ([^(r)]^ + il} + - [Mt)]$_ m+ 

o Jo • 



t A * A 2 

= e_ ^— [ 6 ( T - T s) - S (r - r M J] + e "^ 2(A + A _) " T ^-) - S (t - tm + )] (C54) 

After averaging over tm + and tm_ , the second term will vanish. Since the case Xs < gives the symmetric contri- 
bution via the exchange (A + , A_), we finally obtain that the full correction to the probability distribution reads with 
the measure (|C27|) 

[Mt)]W = [ +C ° — (2@-Y «e-^ X -+ X - [5(t-t s )-5(t-t m _)] » d (C55) 

jR(t) T M \ Tjvf J 2 

We thus finally obtain at first order in /i 

- «(r < *<))£*»> (^) § Kl ( 2 „(0) (M) 4 ) K 2 („(0, 

-o o ,,M0(JiM < t)5 (MY (C56) 

with (JCllf) . The prefactor Qo,i(a*) jQUjl of the second term represents the probability to be in a trap of type S at 
time t and is of order [i l|C32|) . So in the domain t > R, the total distribution [ipt(T)r°' + ^ 1 ' of r keeps the same form 
that at zeroth order, but the amplitude is (1 — f2o,i(A t ))- There is now a contribution of the domain t < R which was 
absent at zeroth order. In particular, in the limit t — > 0, we obtain the essential singularity 

(l) ^ ^4(Q)m e -^mm* ( C57) 

We now compute the correction to the scaling function II(w) (|92|) at first order in \i 

fl (1) («) = H f + dzz 2 K 1 (z)K 2 (z)e' v (^) - Oq.iO) / tfy/V^e"^ (C58) 
So for large v, the only correction to the scaling function IT )(u) l|93|l comes from the second term 



v— »oo 

which corresponds to a correction of order /i 2 for the amplitude of the algebraic decay 1 fv^ 



11™ (v) ~ -Oo.iM-^ (C59) 



b. Probability C(t + t w ,t w ) to be at time (t + t w ) in the same trap at it was at time t w 

Since in our framework, C(t + t w , t w ) is determined by the probability <f> tw (T) to be at time t w in a trap of escape 
time T (|9 f>|) . we obtain as for U(t w + t,t w ) that the configurations M-M+M++ do not give any correction, because 
the escape time of M++ has the same statistical properties as the escape times of M_ and M+ . 

The correction due to the configurations M_5M+ can be written for the case Xs > as 

[MT)]T^ + = e~^ X -~ Xs 5(T Ts) + (1 - e~^) X ~ ~f S 6(T - T M+ ) 



1 — e T s (1 — e T s ) 



X+ + X- 



6{T-T M _) (C60) 
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so that the specific contribution yields after integration over X$ £ [0, A_] 

£~ dX.[Ur)]^_ SM+ = £~ dX. (lMT)]T^ + ~ IMt)]W_m + ) 

= e -^ibz. [ S (T - T S ) - S(T - T M _)] + *~ [S(T - T M _) - S(T - T M+ )} (C61) 

After averaging over Tm + and Tm_ , the second term will vanish. Since the case Xg < gives the symmetric 
contribution via the exchange (A+, A_), we finally obtain that the full correction to the probability distribution reads 
with the measure (|C27ll 



\MT)]M =«/ °° e -^±+±± [d( T-T s )-5(T-T M _)]» d 
Jr. t m \tmJ 2 

n t s \t s J Jo Jo \ + 



R drs_ ( R x " r+ °° ' v 
T s \t~s/ Jo 



+ " dAe - r/Z!L (±Y,l T -, 2R ' T " . I (C62) 



' R T M \ T M , 

so the correction to the scaling function C^h) 197(1 reads 

/■+oo 

&V(h) = I dT{MT)} W e- h ± 
= V 



A 1 \ + +X- 



f *V r d X + r 00 ,A_e-^--)-A_e^ 1+ ^K^ + -) 
Ji v Jo Jo 

-» r ^ r dx+ r ^-^-^-^^ 

+ °°dAe-> f\^e-° h ^<i + ^=) (C63) 



APPENDIX D: SET OF THE IMPORTANT CONFIGURATIONS AT ORDER n 

With the T„ traps described in the text (|102|l . we have to construct the possible f2„ configurations of (2 + rc) traps, 
that are ordered in positions, and that contribute up to order fj, n . We have 

n 

fi n = f2 n _i + uj n = lo % (Dl) 

8=0 

where uj n represents the number of configurations that begin to contribute at order n. We may decompose 

X! 4 Mfc -' Ml+) (D2) 

fc>l,i>l,fe+K2+n 

where the number of configurations that contain Mk- as leftmost trap and Mi + as rightmost trap. For 

(fc = 1, 1 = 1+n), there is only a£ M-,Af(1+n5 = 1 configuration {M_, M + , M2+, ..M(„ +1 ) + }, whereas for (k = 1, 1 = 1), 
there are 

a (M_,M + ) =n , (D3) 

configurations, since we have to order in space the n traps s[°\..S^ in the interval ]M_,M + [. More generally, at 
order (k,l), to construct the configurations of (n + 2) traps containing Mk-M(u_x\_, ..M_, M+, M2+, ■■Mi+j which 
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represent (A; + I) fixed traps, we have to choose (n + 2 — k — I) traps among the (k + I — 1) available intervals and to 
count the possible positional orders in each interval 



+00 +oo /k+l—1 



o£~»- = ^2 ... 6 \ J2 Pi =n + 2-k-l\pi\...pk+i-il 



Pi=0 Pfc + i-i=0 \ i=l / 

The final result is thus that the number of new configurations that appear at order n reads 

>-oo +oo k+l — 1 



(D4) 



E 



k>ld>l,k+l<2+n 



E E § ( E Pi = U+ 2 ~ fc ~ l)pi\--Pk+l-l 

Pl=0 Pfc + !-l=0 i = l 



(D5) 



which generalize what we have found before for the lowest orders with wo = 1 corresponding to (M_, M+), and lq\ = 3 
corresponding to ( M__ , M_ , M+ ) , (M_ , M + ,M++) and (M_ , S, M+) . 



APPENDIX E: USEFUL PROPERTIES OF BESSEL FUNCTIONS 



The following integrals yields Bessel function of type K 



dxx v ' e 



2 [ - ) K 1+ „{2Vab) 



The asymptotic behavior at infinity is independent of v and reads 

K v {z) -< Me~ z 

z-^-oo y ZZ 

Near the origin, the behavior depends on v. We will need the behavior for v = 1 



^0 z 2 



i z 1 

In - + lEuler - ~ 



+ 0{z 3 lnz) 



Another useful integral is 



dzz^Kl{z) = ^ mnk + v)T{k - v) 



AT (k + l) 



(El) 



(E2) 



(E3) 



(E4) 
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